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ABSTRACT 


Single date Landsat CCTs were processed, by Image-100 
to classify Pinus and Eucalyptus species and their age 
groups. The 'study area Mogi-Giiaqu is located in the humid 
sub- tropical climate zone of Sao Paulo State. Prominent 
Pinus species are P. elliottii and P. taeda . Eucalyptus 
species are E. alba and E. saligna . The ages of these 
tree types range fi’om 8 months to 20 years. The study 
area was divided into ten preliminary classes and fea- 
ture selection algorithm was used to calculate 
Bhattacharyya distance (B-distance) between all possible 
pairs of these classes in the four available channels. 

The classes having B-distance values less than 1.30 were 
grouped together resulting in four classes: 1) class PE - 
P. elliottii , 2) class PO - Pinus species other than 
^7 elliottii , 3^ class EY - Eucalyptus spp. under two 
years, and 4) class EO - Eucalyptus spp. more than two 
years. For classification, single cell signature 
acquisition option was used. The low percentages of cor- 
rect classification for class PO was due to its relatively 
small area (a total of 1.44 km2) while for class EO error 
was caused by the similar signature responses of cerrado. 
gallery forest and heterogenous P. elliottii plots in the 
study area. 'The percentages of correct classification 
ranged from 70.9% to 94.12%. Comparisons of acreage esti- 
mated from the Inage-100 with ground truth data showed 
agreemeit. The Image-100 percent recognition values for 
the abovq four classes were 91.62%, 87.80%, 89.89% and 
103.30%. Greater precision of area estimation could have 
been obtained if the ground truth data based on the 
seedling area at planting time had been up-dated. 


1. INTRODUCTION 

Ijandsat is one of the potential tools which is meeting the increasing de- 
mands of government agencies and private industry for timely and precise inven 
torics of forestry data. In the past few years various researchers have tested 
the possibility of identifying forest plant communities using Landsat data by 
conventional or automatic classification techniques. Sayn-V.’it tgensteln- (1972) 
used a color composite of MSS band 4, 6 and 6 to separate coniferous forest 
(mostlji spruce), hardwood (mostly poplar.s) and areas covered by willow a'nd al- 
der. Lopez-Cuervo (1973) visually analyzed Landsat enlarged images and identi- 
fied Pinus and Fiuca l yptus at the ge.neric level according to their tonal differ 
ent iaTions . Gimbarz::vsf:y (1974) also confirmed the separability of coniferous 
and deciduous foi’cst using Landsat imagery. Lee et al. (1974) observed that 
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Landsat imagery could be used to distinguish forested and non-forested lands. 
Difficulties were encountered in separating coniferous species, even though age 
difference was recognized. Kalcnsky and Scherk (1975) reported that overall 
Image-10'0 classification accuracies for coniferous forest, deciduous forest and 
non-forested land ranged from 67% to 81% for single date imagery while for raul- 
tidate imagery consistently above 80%. 

The objective of this study was to classify Pinus and Eucalyptus species 
and their age groups, using the Image-100 automatic classification system. 


2. STUDY AREA AND DATA SOURCES 

Sao Paulo State contributes more than 50% of the paper and cellulose pro- 
duced in Brazil. The state has an artificially reforested area of 6,144 kmZ 
predominantly in pine (22.15%) and eucalypt (76.48%). The study area Mogi-Guaqu 
is located at 22°15' S and 470io' W (Fig. 1) and is a representative of pine and 
eucalypt plantations common to the region. This area includes the Campininha 
pine experimental station of the Forestry Institute of Sao Paulo State (IFSP) 
and the Santa Terezinha eucalypt plantation of the Champion Cellulose & Paper 
Company (CCP). The major Pinus species in Campininha are P. elliottii and 
P. taeda . Other species such as P. caribaea , P. bahamensis , P. oocar pa and 
P. palustr is are also planted in small areas. The prominent Eucalyptus species 
in Santa Terezinha are E. alba and E. saligna . The tree age range from 8 months 
to 20 years. 

Landsat multispectral scanner CCTs of September 13, 1975, orbital number 
178, were processed for the classification of Pinus and Eucalyptus species and 
their age groups by the Image-100 system. 

"Ground truth" information and forest cover maps (Fig. 2) were provided by 
IFSP and CCP and used for supervised classification and accuracy analysis. Spot 
field checks to verify the **ground truth” were also carried out. 


3. IMAGE-100 A.NALYSIS PROCEDURES 


3.1 FEATURE SELECTION 

Prior to the supervised signature extraction and classification, the CCTs 
were corrected radiometrically . A preliminary analysis of the study area was 
carried out through four B & W imagery and a false color composite displayed on 
the image monitor at an approximate scale of 1:82,000. The following ten pre- 
liminary classes were dei'^ived for the Image *l00 study. 

PE_ - P. elliottii over 13 years. 

d - " ' 

PEjj - P. elliottii under 13 years. 

OP - Pinus species other than P. elliottii and P. taeda . 

PT - P. taeda . 

^8mo 2yrs ” Eucalyptus spp. from 8 months to 2 years. 

^2yrs. ~4yrs. “ Eucalyptus spp. from 2 to 4 years. 

^4yrs. - 7yrs. " spp. from 4 to 7 years. 

^7yrs. - 20yrs. ' Eucalyptus spp. from 7 to 20 years. 

RE ^ 2yrs ~ Eucalyptus spp. under 2 years. 

^^2yrs - Gyr*" Eucalyptus spp. from 2 to 9 years. 


Bhattacharyya distance (B-di stance, Mari 11 and Green, 1963) between these 
ten preliminary classes war, calculatfd using INPF.’s feature selection algoritlim. 
This algorithm assumcr: that each class has a multivariate Gaussian distribu- 
tion. In this operation, 40 to 50 tralTvipg samples of each class were deline- 



ated on the image monitor by the electronic cursor which was adjusted to the 
minimum size equivalent to a pixels/ sample. Then the B-distance value was 
calculated in the four available channels. The results of the B-distance values 
for the preliminary Pinus and Eucalyptus classed in four channels are presented 
in Tables I and II. The separability of the preliminary classes was defined 
according to the values of B-distance. Using the curve of Swain and King (1973) 
the estimated probability of correct classification is > 0.8S for the value of 
B-distance > 1.30. Thus those classes having B-distance values less than 1.30 
were grouped together to form four classes as follows: 

PE - — P. elliottii . 

PO --- Pinus species other than P. Elliottii . 

EY --- Eucalyptus spp. under two years. 

EO — - Eucalyptus spp. over two years. 

Table III lists the "B-distance" values of classes PE» P0» EY and EO in 
four spectral channels. 

3.2 SIGNATURE EXTRACTION AND CLASSIFICATION 

Single cell signature acquisition of the Image-100 was used for classifi- 
cation. This option creates a four-dimensional parallelepiped; each of the sides 
of which corresponds to the signature limits of the training areas in each 
channel. Training samples, independent from the test area, were selected from 
each of the four classes (PE, PO, EY and EO) using the electronic cursor. The 
limits of the histograms of the training areas were modified in the four 
spectral channels until a satisfactory classification was obtained. The spec- 
tral statistics of the histograms .from the training samples of four classes are 
given in Tables IV to VII, respectively. A density slicing routine was also 
used to give pixel frequencies of the study area for all possible combinations 
of the four channels selected axes. An example of this result is shown in Fig. 3. 


4. IMAGE-100 PERFORMANCE 
.1 PIXEL CLASSIFICATION ACCURACY 

After classification, the electronic cursor was again positioned on the 
image monitor. The "alphanumeric theme print" was executed for a cursored area 
of approximately 50 km2. This printout was compared pixel by pixel with the IFSP 
and CCPs forest cover maps to produce a confusion matrix (Table VIII). This 
table shows the number of pixels classified according to its "true" class and 
Image-100 *s "chosen" class in the cursored area. The percentages given in paren- 
thesis were calculated based on the total number of classified pixels in each 
class. Histogram modification thresholding in the four spectral channels tends 
to increase errors of omission and decrease errors of commission. Thus, percent- 
ages of correct classification which take into account errors of omission were 
not high. To give a more realistic view of the Image-100* s pixel classification 
accuracy the following formulas were used to calculate the percentages of cor- 
rect classification and errors of commission. 


Correct Classification (%) 


No. of correctly classified pixels in class A 
Total no. of pixels in the test area of class A 


Error of Commission (%) 


No. of pixels classified as class A but actually 

belonging to another cl.iss 

Total no. of pixels in the test area of non-A class 


4.2 AREA ESTIMATION ACCURACY 

The Image-100 estimated class areas were compared to their corresponding 
"ground truth" data provided by IFSP and CCP. The pcrceiitagc of the estimated 
areas on ground truth data are presented in Table X. 
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5. DISCUSSION AND CONCLUSIONS 


For the study, sufficiently large test areas of > pine trees less than 
thirteen yeai^’s old and trees over thii’teen years old were further separated 
geographically into two classes, PE^ and PE^,, although a ”B-distance" value of 
0.22 indicated that no satisfactory separation was possible. P. taeda and the 
class OP were readily separable from P. elliottii but unseparable from each 
other (Table I). The morphological similarity among the Eucalyptus species 
made their differentiation difficult, although Image-100 successfully classi- 
fied Eucalypt us under two years from the other preliminary age groups (Table II). 
This separabiTity was due to the distinct spectral characteristics of the homo- 
genous young Eucalyptus which aided the Image-100 identification. 

For classification a hyperparallolepiped signature file was constructed 
for each of the classes PE, PO, EY and EO using a "single cell signature acqui- 
sition" program. The spectral statistics of the signature files are shown in 
Tables V, VI, VII and VIII, but separability of these four classes was most 
easily observed on "N-dimensional histogram two-dimensional print" (Fig. 3). The 
longer distance between genera Pinus and Eucalyptus confirmed their 
highei' separability than species and age classes. The percentages of Image 100 's 
correct classification for the classes PE, PO, EY and EO were found to be 

89.14%, 79.26%, 94.12% and 70.99% respectively. These percentages were rela- 

tively lower compared to those calculated using the total number of pixels from 
each class classified as the denominator (Table IX). The lower percentage of 
correct classification was obtained from the class EO (70.99%). It was caused 

by cerrado, gallery forest and heterogeneous P. elliottii plots of the area 

which gave a similar spectral responses as class EO. The highest percentage of 
correct classification was found for class EY (94.12%). It may have been due to 
the homogeneity of the eucalypt under two years old. 

In the area estimation study the misclassified pixels shown in Table VIII 
accounted for the low pixel classification accuracies of Table IX and gave area 
estimations of 17.28, 1.44, 4.00 and 22.51 km2 for the classes PE, PO, EY and 
Eo respectively (Table X). The -above areas estimated by the Image-100 when com- 
pared to the same ground truth areas gave recognition values (RVs) of 91.62%, 
87.80%, 89.39% and 103.30%. The 3.30% overestination of class EO was due to a 
high proportion of misclassified pixels (13. 61%) v;hich had produced spectral 
signature similar to class EO. The overall area-weight difference of 6.31 was 
achieved. The results of the Image-100 area estimations were encouraging, con- 
sidering that the machine classification was carried out at the species level, 
for pine, and based on age groups, for eucalypt. 

The ground truth provided by IFSP and CCP agenci es were based on the 
seedling area at planting time. These data represented "theoretically refot-ested 
areas" under "optimum conditions" without damage and mortality. In the field 
check, various plant-failure areas were observed. The recognition value could 
have been greater if the ground truth data had been up-dated annually. 

The following significant results were obtained from automatic classifica- 
tion of reforested areas of pine and eucalypt: 

Based on "B-distance" values the following spectral separations were achieved: 

1, P. elliottii was separated from the ot; ir pine species. 

2* Eucalyptus spp. under two years old were separated from 
Eucalyptus spp. more than two years old. 

3. Clear distinctions v;ere made between pine and eucalypt. 

Based on Image-100 area estimations, RV ranged from 87.80% to. 
103.30% for the four classes defined (PE, PO, EY and EO) when compared to the 
same information provided by forestry agencies. 

These results point to the operational potential of computer automated 
processing of Landsat data to recognise tree species and calculate standing 
volume, especially in remote regions whei-e gi-oui.d and aerial surveys arc often 
costly. 



Before the methodology used in this study can be applied on an operational 
scale, more consideration should be given to signature extention to other geo- 
graphical regions, and to spatial differences between areas of natural and ar- 
tifical reforestations. 
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Figure 1. LOCATION OF STUDY AREA IN MOGI-GUAQU, SAO PAULO STATE 
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Figure 3. 


AN EXAMPLE OF ’’N-DIMENSIONAL HISTOGRAM TWO-DIMENSIONAL PRINT" 
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Figure 5. ALPHANUMERIC THEME PRINT USED FOR PIXEL CLASSIFICATION AC- 
CURACY MEASUREMENT 
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TABLE IV. SPECTRAL STATISTICS FROM THE TRAINING SAMPLES OF CLASS PE 


CHANNEL 


SPECTRAL 

BOUNDS 

LB 

UB 


DELTA 

MEAN 

2 

18. S 

S 

13.3 

£ 

19.1 

3 

25.4 


VARIANCE 


.2 


TABLE V. SPECTRAL STAT'STICS FROM TKI TRAINING SAMPLES OF CLASS PO 


CHANNEL 


SPECTRAL 

BOUNDS 

LB 

UB 

18 

21 

12 

19 

35 

42 

42 

56 


DELTA 


VARIANCE 


TABLE VI. SPECTRAL STATISTICS FROM THE TRAINING SAMPLES OF CLASS EY 
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TAbuE VIII. CONFUSION MATRIX 
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TABLE X. COMPARISON OF IMAGE-100 ESTIMATED AREAS WITH "GROUND TRUTH” DATA 










































